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■ We present a simple way to quantize the well-known Margulis expander map. The 
^ l' result is a quantum expander which acts on discrete Wigner functions in the same way 
^ I ^ ^ the classical Margulis expander acts on probability distributions. The quantum version 
^ , shares all essential properties of the classical counterpart, e.g., it has the same degree and 

■ spectrum. Unlike previous constructions of quantum expanders, our method does not 
'~{ ' rely on non-Abelian harmonic analysis. Analogues for continuous variable systems are 

mentioned. Indeed, the construction seems one of the few instances where applications 
I I , based on discrete and continuous phase space methods can be developed in complete 

1 analogy. 
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\^ , Motivated by the prominent role expander graphs play in theoretical computer science 

' [1], quantum expanders have recently received a great deal of attention [21 [3l HI [3 [H [7]. 

■ In this work, we report an observation which allows for the simple explicit construction of 
such quantum expanders. The method relies heavily on quantum phase space techniques: 
Once familiar with this techniques, the result is an almost trivial corollary of the analogous 
classical statement. We further discuss continuous analogues of quantum expanders, where 

. 5^ I again, phase space methods render this an obvious generalization. Hence, the present work 

' can equally be regarded as the presentation of a simple quantum expander, as as a short 

^ , exposition of the strengths of the phase space formalism as such. 

1 Preliminaries 
1 . 1 Expanders 

Expander graphs turn up in various areas of combinatorics and computer science (for all 
claims made in this section, the reader is referred to the excellent survey article Ref. jl]). 
They often come into play when one is concerned with a property which "typically" holds, 
but defies systematic understanding. A simple example is given by classical error correction 
codes. One can show that a randomly chosen code is extremely likely to have favorable 
properties, but it seems very difficult to come up with a deterministic construction of codes 
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Fig. 1. The phase space distributions resulting from three appHcations of the Margulis expander 
acting on a configuration initially concentrated at the origin of a 7 X 7 lattice. The starting 
distribution can be interpreted either as a classical particle with a well-defined position on a 
two-dimensional lattice, or as the quantum phase space operator A(0, 0) (see text for definition). 

which are "as good as random" . Expander graphs can be expHcitly constructed, but capture 
some aspects of generic graphs. It turns out that this property can be used to de-randomize, 
e.g., the construction of codes or certain probabiHstic algorithms. 

The formal definition is straightforward. Consider a graph G with N vertices V , each 
having D neighbors (we allow for multiple links and self- links). There is an obvious way to 
define a random walk on the graph: At each time step, a particle initially located on a vertex 
V will be moved to one of the D neighbors of v with equal probability. The resulting Markov 
process is described by an x iV doubly stochastic matrix A. The largest eigenvalue of A is 
Ai — 1, corresponding to the "totally mixed" eigenvector (1, . . . , 1)/N. Let A be the absolute 
value of the second largest (by absolute value) eigenvalue. A small value of A means that the 
Markov process is strongly mixing, i.e., converges rapidly to the totally mixed state. We call 
G an (A^, D, A) expander if it is described by these parameters. The goal is to find families of 
expander graphs with arbitrarily many vertices A^, but constant (and small) degree D and A. 

While the notion of an expander graph seems hard to quantize (see, however, Ref. |3]), 
it makes sense to look for quantum analogues of strongly mixing Markov processes with low 
degree. Indeed, we call a completely positive map A a (N, D, X)-quantum expander if A can 
be expressed in terms of D Kraus operators acting on S((D^) and the absolute value of its 
second largest singular value is bounded from above by A (here, B{T-C} denotes the space of 
linear operators acting on a linear space 7i). Once more: The intuition is to have a quantum 
channel which can be written using few Kraus operators, but which rapidly sends any input 
to the completely mixed state under repeated invocation [8]. 

Quantum expanders have been introduced independently in Ref. [2] in the context of 
Hamiltonian complexity for the purpose of constructing states of spin-chains with certain 
extremal entanglement and correlation properties, and in Ref. [3], where the problem was 
approached from a computer science perspective. Very recently, randomized [3] and explicit 
[H HI El IS] constructions of expanders have appeared in the literature. The basic idea of a 
quantum expander is implicit in earlier work [7j. 

1.2 Margulis expander 

Margulis provided the first explicit construction of a family of expander graphs [S]. Their 
expansion properties can be verified by elementary (if tedious) means [1]. 
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The vertices of Marguhs' graph are given by the points of a iV x A^-lattice|"We label the 
axes of the lattice by the elements of TLm — {0, . . . , iV — 1}. Now consider the four afline 
transformations on given by 

Ti : v^Siv, (1) 

T2 : Siv + {l,Qf, 

T3 : V ^ S2V, 

T4 : v^S2V-[Q,lf, 

where 
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All operations are modulo N . Let S be the set of these four operations, together with their 
inverses. In Margulis' construction, two vertices are considered adjacent if and only if they 
can be mapped onto each other by an operation in S. 

One finds that A is bounded above by \/2 5/8, independent of N [101 [T]. An instance of a 
random walk on the Margulis graph is visualized in Fig. [1] 

1.3 Discrete phase space methods 

The present section can be approached from a purely mathematical, or from a physically- 
oriented point of view. To make the argument more accessible, we will briefly outline both 
approaches before going into details. 

Mathematically^ one starts by noting that the operations used in the construction of 
the Margulis graph (Eq. ^) are afflne transformations on Z^, where in addition the linear 
part Si has unit-determinant. The functions of this kind form a finite group, which we will 
refer to as Gm- Two facts will be established below. Firstly, 

• there is a (projective) unitary representation 

T^Ut 

of Gn on 

This representation facilitates the quantization of the expander. Indeed, the quantum Mar- 
gulis expander will be defined as the c.p. map A^r, which applies one of the unitaries [/t, T € S 
at random. To prove that this construction defines an expander, we will need a second fact: 

• Let N be odd. There are iV^ hermitian operators A{a) G B{(D^), labeled by the points 
a G Z^, such that 

1. The operators form an orthonormal basis with respect to the Hilbert-Schmidt inner 
product: 

^ti-{A{v)A{w)) ^6^,^. (2) 



" Note a slight inconsistency in our notation: the number of vertices in this case is N'^, not A'^. 
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2. The basis thus defined is compatible with the unitary representation of Gjq in that 

UTA{v)u}r^A{T{v)), (3) 

for T e Gat, v G TL%. 

In order to analyze the action of Kn on a density operator p, we will use Eq. ([2|) to expand 
p in terms of the A(a)'s and then Eq. ^ to reduce the problem to the classical case (see 
Section El). 

Physically speaking, we will employ a phase space description of the quantum system. 
Recall that the term phase space originates in classical mechanics. Here, the state of a 
single particle in one spatial dimension is completely specified by two real parameters: its 
position and its momentum. The two-dimensional real vector space spanned by the position 
and the momentum axes is referred to as the particle's phase space. Likewise, the state 
of a single quantum system can be specified by a quasi-probability distribution on phase 
space, namely the particle's Wigner function. The Wigner function shares many properties 
of classical probability distributions, except for the fact that it can take negative values (see 
Ref. [TTJ[T21[I3] for an analysis of quantum states which exhibit only positive values). 

In the context of continuous-variable systems, affine volume-preserving transformations 
of the phase space are known as canonical transformations. Let p be a density matrix and 
denote by Wp{v) the associated Wigner function. It is well-known [5T] that for every canonical 
transformation T, there is a unitary operator Ut which implements the map T in the sense 
that 

As detailed below, a similar relation holds for finite-dimensional quantum systems, associated 
with discrete phase spaces. Indeed, the Wigner function of a density operator p turns out to 
be nothing but the collection of expansion coefficients of p with respect to the basis given in 
Eq. ((2|); canonical transformations are elements of Gat; and the correspondence T Ut is 
just the representation mentioned in the first paragraph of this section. 

So in this physical language, the basic realization is that the building blocks of the Margulis 
scheme (Eq. ([l}) are canonical transformations of a discrete phase space. 

To make all this more precise, let iV be od(^ Ti- — (D^ and assume that some basis 
{|0), . . . , |iV — 1)} in Ti has been chosen. Let lo = eT?* be an A^th root of unity. We define 
the shift and boost operators as the generalizations of the X and Z Pauli matrices by 

x{q)\k) =\k + q), z{p)\k) ^ uP'^\k) (4) 

(arithmetic is modulo N). The Weyl operators are 

w{p,q)^LO-^''P^z{p)x{q), (5) 

where = {N + l)/2 is the multiplicative inverse of 2 modulo N . For vectors a = (p, q) G 
7L\i, we write w{a) for w{p,q). Let 

A{0,0) -.{x) ^\~x) (6) 

' We restrict attention to odd dimensions, as the tlieory of discrete Wigner functions is much more well-behaved 
in this case. 
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be the parity operator and denote by Aija^ q) its translated version, 

A{p,q)^w{p,q)A{0,{))w{p,q)^. (7) 

We will refer to the A(p^qys as phase space operators. One can check by direct calculation 
that Eq. ^ holds. The Wigner function of an operator p is the collection of the expansion 
coefficients of p with respect to the basis formed by the phase space operators. Formally: 

W,{p,q) = ^tr{A{p,q)p). (8) 

There are two symmetries associated with a phase space: translations and volume-preserving 
linear operations. We shortly look at both in turn. Firstly, it is simple to verify that for 

w{a)A{b)w{a)''^A{a + b). (9) 

Hence, Weyl operators implement translations on phase space. Secondly, let 5 be a unit- 
determinant matrix with entries in W^n- It turns out [171 1151 [T^ that there exists a unitary 
operator /i(S') such that, for all a G the relation 

p{S) A{a) ^i{Sy ^ A{S a) (10) 

holdflf 

It follows immediately that for every affine transformation T of the type given in Eq. ([T]) , 
there exists a unitary operator Ut such that 

WurpU^i^) = Wp{T~\a)). (11) 

Hence, one can unitarily implement the building blocks of Margulis' random walk. 

2 A quantum Margulis expander 

With these preparations, it is obvious how to proceed. Define the completely positive map 
An by 

AN{p)^^^Y.UTpUi., (12) 
' ' Tes 

where we have used the notation defined in Eq. pT|) above. One immediately gets: 

Observation 1 (Quantum Margulis expander) For odd N, the map (Eq. hl^) ) acts 

on Wigner functions in the same way the Margulis expander acts on classical probability 

distributions. In particular, its degree and its spectrum are identical to the ones of the Margulis 

random walk. The Wigner functions of A's eigen- operators are the eigen- distributions of the 

classical random walk. 

(c) 

Proof. Let Aj^^ be the stochastic matrix associated with the random walk on the classical 
Margulis graph. For v e K^, let e{v) be the function on K^, which takes the value 1 at w 
and else. Clearly, the set {e(w)}t,g^2^ spans the space of all functions on the lattice. Also, 

A'^Pieiv)) = ^\2<^(^)y 

The operator /i(S') is the metaplectic representation of the symplectic matrix S. In quantum information 
theory, the set {w{a) fJ,{S) : a £ Z5^,det(5) = 1} is called the Clifford group 1141 . which must to be confused 
with the Clifford group appearing in the context of Fermions or representation theory of SO{n). 



6 Quantum Margulis expanders 

Using Eqs. (fTl] [T2 |) . we get for the quantum version 

Aa'(A(z;)) = ^^A™). 
' ' Tes 

Hence the action of the classical and the quantum expanders are identical on a basis 
claims follow. 

3 Efficient implementation 

Consider a quantum expander which acts on a tensor-product Hilbert space ((C*)®" ~ (D^ 
for N — d". The expander is efficient if it can be realized using poly(rt) single-qudit or 
two-qudit quantum gates. So far, only two efficient constructions have been published [51[B]. 
The Margulis expander adds to this list. 

Tfieorem 2 (Efficient implementation) The quantized Margulis expander acts effi- 
ciently on (C'')®". 

To establish the claim, we need to clarify how we introduce a tensor product structure 
in <D^ . Every < j < — 1 can be expressed in a d-adic expansion as j = ji . . . jn for 
< j; < d. More precisely, j = 'Ym=i ■ The tensor product structure is now given by 

\j) = |ji)(g)---(g)|j„). 

Lemma 1 (Efficient constituents) Let N — d". The following operators act efficiently on 

1. The quantum Fourier transform 

F:|j)K^iV-i/2^exp(z— j-A;)|fc). 

fc=0 

2. The Weyl operators w{l,0) andw{0,l). 

3. The operators fJ-iTi) and fi{T2). 

Proof. The first statement is well-known. See Chapter 5 in Ref. [18] for the qubit 
version, which can easily be adapted to general d. Next, consider w{l, 0) — z{l). We have 

z{l)\j) = exp(i^j)|ji,...,j„) 

71 

= exp(i27r ^j;d"')|ji,...,j„) 

1=1 

(g)exp(i27rj,d-')bi). 

Hence is actually local. One confirms that x{l) = F z{l) F'' and thus a;(l) is efficient. 

To conclude the proof, we need to borrow three statements from the theory of metaplectic 
representations. Firstly, is a projective representatior|f i.e., 

^l{ST) = e"^(^'^) t^{S)n{T) 

''Actually, fj, is even a faithful representation, but that fact is irrelevant for our purposes. 
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for some phase (j){S,T) [H [HI US] . Secondly, 



1 
-1 



and thirdly, 



is given by 



u± = ^i{ 



1 ±2 
1 



f/±|j) =exp(i2^/iV(T/))|j). 



The last two statements can be found in Theorem 4.1 of Ref. 15J (strictly speaking only for 
the case of prime iV, but the proofs work for any odd value) or in Lemma 2 to Lemma 4 of 
Ref. [17]. The claim becomes easy to verify: 

n 

U±\j) = exp(z2^(T 5] jzjz'd"'''''))b> 

= \{R{l,l')\3), 
1,1' 

where we have introduced the diagonal two-qudit unitary 

R{l,n\ji,Ji') = exp(z2^(TjzJi'd"-'"'')) \ji,Ji')- 
Thus U± - and therefore in particular /i(Ti) - are efficient. Finally, 



T2 = 



" 1 ■ 




■ 1 -2 ■ 




■ 1 " 


-1 




1 




-1 



which implies that /i(r2) oc FU-F ^ is efficient. q 

The proof of Theorem[2]is now immediate, as all the ?7t's which appear in the construction 
of A can be implemented by combining the unitaries treated in the above lemma and their 
inverses. 



4 Continuous variable systems 

The quantum phase space terminology of Section 11.31 has originally been introduced in the 
context of continuous variable systems (see e.g. Ref. Ji) . In particular, if we re-interpret 
the affine transformations S given in Eq. ([T]) as operations on IR^, we immediately obtain a 
completely positive map Aqo acting on the infinite-dimensional Hilbert space of a single mode. 
Does it constitute a quantum expander? After reviewing some definitions in Section [4.11 we 
will give an affirmative answer in Section [4.2l The action of expanders on second moments is 
discussed in Section [473l 
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4-1 Continuous phase space methods 

In the continuous case, the phase space is given by R^. Let X and P be the canonical position 
and momentum operators. The Weyl operators [191 UHl HI] are now 

w{p, q) — exp(i(7P — ipX). (13) 

As in Eq. the parity operator A{0, 0) acts on state vectors ijj E L^(]R) as 

{A{0,0)i;){x) ^ i;i-x). 

We define the phase space operators A{p, q) for (p, q) e exactly as in Eq. ([7]). The Wigner 
function becomes 

Wp{p, q)^Tr-Hr (A(p, q) p) 

c.f. Eq. (O. The obvious equivalents of Eqs. (jOTTOl) hold for a G R^ and 5 g Sp(2,R), 
the group of unit-determinant transformations of the two-dimensional real plane. Hence it is 
plain how to interpret Eq. pT|) and finally how to turn Eq. JTS]) into a definition of Aqo, the 
infinite-dimensional quantum Margulis map. 

4.2 A continuous quantum Margulis expander 

A slight technical problem arises when transferring the definition of an expander to the infinite- 
dimensional case: both the invariant distribution f(v) = 1 of a classical expander and the 
invariant operator 1 of a quantum expander map are not normalizable. Hence, if we define 
e.g. the action of a completely positive map A on the set of trace-class operators T^{'H), 
the would-be eigenvector with eigenvalue 1 is not even in the domain of definition. In the 
light of this problem, we switch to the following definition of a quantum expander, which is 
compatible with the notion used up to now. 

Definition 1 Let N < 00 and set Ti, = (D^. A completely positive map A is an {N,D,X)- 
quantum expander if, for all traceless operators X G T^(7i), 

||A(A)||2<A||A||2. 

The definition above is best understood in terms of the Heisenberg picture: 

|tr(A"(p)A)| = |tr(p(At)«(A))|<A" 

for all normalized (| 1^^112 = 1), traceless observables X. Thus the state becomes "featureless" 
exponentially fast when being acted on by A. Let Xm be the second largest eigenvalue of the 
finite Margulis expanders. Then: 

Observation 3 (Continuous quantum expander) The infinite- dimensional quantum Mar- 
gulis map Aoo is an {co, 8, Xm)- quantum expander. 

Note that by the previous section, we know there are {N, 8, Xm) quantum expanders for 
arbitrarily high N. A priori, however, this does not imply the existence of a solution for 
A^ = 00. 

Once more, by switching to the phase-space picture, the proof of Observation [3] can be 
formulated completely in classical terms. The intuition behind the argument is simple to 
state. Take an element T of S, e.g. 

T:v^\l l]v. (14) 



D. Gross and J. Eisert 9 



The inverse is given by 

' 1 -2 
1 



(15) 



regardless of whether the matrix is interpreted as acting on ]R^,Z^ or Z^. As the same 
is true for all other elements of <S, the action of the classical Margulis map "looks similar" 
on continuous, infinite discrete and on finite phase spaces - at least as long as it acts on 

distributions which are concentrated close to the origin, so that the cychc boundary conditions 
of do not come into play. Using this insight, the following lemma reduces the continuous 
to the finite case. 

Lemma 2 Let f G Co(]R^) be a continuous function with compact support, such that 

f{v)dv = 0. (16) 
Let A : i^(R^) — > L^(]R^) he the classical Margulis map acting on distributions on R^. Then 

||^(/)||2<Am||/||2. (17) 

Proof. We discretize the problem by partitioning into a net of squares with side 
length 5. More specifically, for {x,y) G 1?, let 

Q,{x, v) = [{x - 1/2)5, {x + 1/2)5] x [{y - 1/2)5, [y + 1/2)5] 

be the square with edge length 5 centered around {x5,y5) G 'E? . The discretized version of 
/ is /a : ^2 ^ (D defined by 

fs{x,y) = 4 / f{v)dv. 
" jQs(x,y) 

Note that Ylix,y fsi^^ v) = 0- On we use the ^-dependent norm 

\\m2 = {5'Y.\f^{x,y)\^f" 

(the factor 6"^ corresponds, of course, to the volume of the squares Qg{x,y)). Now, let T be 
one of the affine transformations in <S. We can interpret T as an operation on Wp and define 
its action on fg accordingly by 

{T{fs)){x,y) = fsiT-\x,y)). 

For small enough 5, the approximation is going to be arbitrarily good: using the uniform 
continuity of /, and the fact that all T E S are continuous and volume-preserving, one finds 
that for every e > 0, there is a 5 > such that simultaneously 

IIIMI2-II/II2I < e/2, (18) 
|||^(/,)||2-||A(/)||2| < e/2. (19) 

As the support of / is compact, there is an ii e N such that fs{x,y) and A{fs){x,y) are 
equal to zero whenever |x| > R ov \y\ > R. This enables us to pass from to the finite 
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lattice Wiff for N > 2R. Indeed, when we re-interpret fs as a function — > (D and the T e 5 
as affine transformations on Z^, the values of ||/5||2 and ||^(/5)||2 remain unchanged. But 
we know that A is an {N, 8, AA/)-expander for every finite N. Hence 

11^(^)112 <Am||/^||2, 

implying (by Eqs. (|18ll9p ) 



II^(/)I|2<Am| 

This proves the claim, as the right hand side can be chosen to be arbitrarily small. |— | 

Proof (of Observation Once again, the quantum Margulis map Aqo acts on the 
Wigner function Wx of any operator X in the same way the classical Margulis scheme acts 
on distributions on R^. Now, X e T^{n) imphes Wx & L^(IR^). Because C^(IR^) is dense 
in L^(]R^) and Aoo is continuous, Lemma [2] suffices to establish the claim. q 

4-3 Action on second moments 

In physics, one often measures the concentration of a phase space distribution by its second 
moments with respect to canonical coordinates. Thus, it may be interesting to look for 
signatures of the strong mixing properties of a quantum expander in its action on second 
moments. 

More precisely, first moments are the expectation values of the position and momentum 
operators {{X),{P))'^ (where (A) = ti{pA) for an operator A). The second moments are 
defined as the entries of the covariance matrix: 



7 = 2Re 



{X^)-{X)^ {XP)-{X){P) 
{PX)-{X){P) {P')-{P? 



As the action of the continuous quantum expander in state space is defined via the metaplectic 
representation, the change in second moments can be computed explicitly. In particular, any 
S £ Sp(2,]R) gives rise to a congruence 7 1-^ S'yS'^ for second moments. More generally, 
it is not difficult to see that for arbitrary convex combinations of states subject to affine 
transformations, the output's first and second moments depend only on the same moments 
of the input. 

Under the Margulis random walk, one obtains for the first moments 

with {xt,Pt)^ — T{{X), (P))'^ ■ For the second moments: 

2 

7 ^ fil) E (^'^^^ + T-'j{T-T) + 2 G, (20) 



where the matrix G is given by 
G 



T |5| VzL^T |5| |5| /-^T,T' |5|2 ' 

T |5| l^T,T' |5P \S\ \l^T \S\) 
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The latter matrix is evidently positive [22j. To show that the main diagonal entries of /^"■'(7) 
diverge exponentially in the number n of applications of the map /, it is hence sufficient to 
consider the map 

7-5(7)-E(t.7t;^+7;-S(t;-T), 



i=l 



Since 



/(")(^)>gW(^) 



A simple calculation yields 



a b 
b c 



5(7) 



a + 2c b 
b c + 2a 



Let 7^") = g'^^\l) be the covariance matrix after n iterations of g and define a ^ [a + c)/2, 
and (3 = (a — c)/2 to simplify notation. Then 



v(") - 



3"a + (-l)"73 b 

b 3"a-(-l)"/3 



This means that 



-logs (7^"^) 



1 
1 



(n — > 00). 



Thus, the elements of the main diagonal - and therefore also tr(/'^"^ (7)), det(/'^"^ (7)), and 
spec(/^"-' (7)) - diverge exponentially in the number n of iterations. 

5 Summary and Outlook 

Employing phase space methods, we were able to quantize a well-established combinatorial 
structure with almost no technical effort. Until now, discrete Wigner functions have been 
studied mainly for their mathematical appeal. As far as we know, the present work is the 
first instance where a problem not related to the phase space formalism itself has been solved 
using the properties of discrete Wigner functions. 

The unitaries which appear in the construction of expanders have randomization properties 
which are in some sense extremal. It would be interesting to see whether connections to 
other extremal sets of unitaries - e.g., unitary designs |23[ I24j - can be found. Also, more 
practical applications may be anticipated, e.g., when one aims at initializing quantum systems 
in the maximally mixed state with few (i.e. D) operations, under repeated invocation of 
the same completely positive map A. From a physics perspective, the existence of quantum 
expanders gives rise to the quite counterintuitive observation that small environments can lead 
to strongly mixing and decohering dynamics: One may think of a controlled-[/ operation, the 
unitaries being those from a quantum expander, where the control has the dimension of a mere 
log2(-D). Lastly, the programme may improve the understanding of iterated randomization 
procedures, as the one discussed in Ref. [25^. 
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